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In this paper, an output-feedback-type wide area damping control (WADC) strategy is presented as a 
supplementary control of high-voltage direct current (HVDC) transmission to enhance thestability of 
large-scale interconnected power systems. This strategy can effectively suppress the varying-delay 
effects of wide area signal on the damping performance. A discrete-time system model that considers 
the varying-delay feature of wide area signal is established, then a dynamic output-feedback (DOF) type 
WADC is formulated by employing a new stabilization analysis method based on Lyapunov function. All 
the solutions of controller design problems are converted into the framework of linear matrix 
inequalities (LMIs). Compared to the conventional design method, the varying-delay of wide area 
signal is involved during each step of controller design, which is advantage to eliminate the effects of 
signal varying-delay on damping performance. Both the four-machine two-area test system and the 
sixteen-machine five-area interconnected system are used to validate the control concept and the 
designed controller. 

© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Nowadays more and more high-voltage direct current (HVDC) 
transmissions are being applied in power systems for long- 
distance power transmission and interconnection of electric net- 
works. As a backbone interconnected line, HVDC transmission 
plays an effective role on system operating performance. It can 
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Nomenclature 


CCL cone complementary linearization 
DOF dynamic output-feedback 

FACTS flexible AC transmission system 
HLF high-pass filter 

HVDC high-voltage direct current 

LFO low-frequency oscillations 


LMIs linear matrix inequalities 

NETS New England Test System 

NYPS New York Power System 

WAC wide area control 

WADC_ wide area damping control 
WAMS wide area measurement system 
PSS power system stabilizer 

SS state-space 


SOF static output-feedback 

TF transfer-function 

x(k) state vector 

u(k) control-input 

y(k) control-output 

A system matrix 

B input matrix 

Ca output matrix with delay effect 


d(k) time-varying delay 

AcBoCoDe state-space matrices of DOF controller 
Ksor control gain of SOF controller 

PLL real matrix 

Q; (i=1,2,3) real matrix 

ZR; G=1,2) real matrix 

X,Y,N,S,M real matrix 


contribute transmission capacity increase, power flow regulation, 
and stability enhancement for interconnected power system 
[1-3]. Meanwhile, with the technical development of synchro- 
nized phasor measurement, wide area measurement system 
(WAMS) gets more and more application in large-scale power 
system [4,5]. WAMS can not only realize global state estimation 
for power system, but also be construct wide area control (WAC) 
strategy by introducing wide area signals [6]. 

In large-scale interconnected systems, such as the European 
interconnected networks [7], America interconnected systems [4], 
and power grids of China [8], inter-area oscillations among 
different areas are easily excited by disturbances, faults or changes 
of system operating condition. This kind of low-frequency oscilla- 
tions (LFO) generally represent the interaction among various 
generators located in different areas, and further represent the 
power oscillations on backbone interconnected lines, thus inevi- 
tably endanger the stability and security of system. 

In recent years, several control methods were proposed to 
prevent LFO, such as the wide area power system stabilizer (PSS) 
[9] located at the generation side, and the flexible AC transmis- 
sion system (FACTS) devices [10,11] placed at the transmission or 
the distribution sides. In practice, HVDC transmission is generally 
used as a backbone line in large-scale interconnected systems. It 
is able to naturally enhance system stability. Further, if the 
supplementary control function of HVDC transmission is devel- 
oped to construct oscillation damping strategy where the control- 
input is provided by WAMS, then the so-called HVDC wide area 
damping control (WADC) strategy could contribute to the effec- 
tive LFO damping effectively and thus improve the overall 
stability of interconnected system. 

However, the delay effect of wide area signals should be 
considered carefully for the WAC strategy. During the period of 
signal transmission and processing, it inevitably leads to a certain 
delay for wide area signals. Due to the technical development of 
WAMS, the time delay can be controlled in a relatively small value 
such as 300ms, but it still affects even destroys the control 
performance [9,12], especially on the practical situation the delay 
is not constant but time-varying. To overcome this problem, some 
control methods that consider the delay effect were proposed [9,13]. 
In fact, power system with the application of WAMS can be looked 
as a delay-dependent system defined in the field of control theory 
[14,15]. It is feasible to employ advanced control theories or 
methods to solve the delay effects of wide area signals, and design 
an advanced controller to fundamentally suppress delay effects. 

In this paper, a dynamic output-feedback (DOF) WADC for 
HVDC transmission will be proposed. The hybrid HVDC/AC 


interconnected system will be linearized as a discrete-time model 
with the varying-delay feature of wide area signal. Further the 
control problem will be formulated as a solution of linear matrix 
inequalities (LMIs). Both the four-machine two-area test system 
and the sixteen-machine five-area interconnected system will be 
used to validate this control design method. 


2. Description of HVDC wide area damping control 


Fig. 1 shows the basic structure of HVDC WADC strategy. The 
rectifier side has a local pole controller and a wide area controller 
which are used to implement the constant DC current control and 
the oscillation damping control, respectively. The WADC is in fact 
a supplementary control of HVDC local pole control, which has 
some features as follows: 


e To get a good damping performance, the WADC chooses 
remote variable (e.g., current or power flow on the tie-line) 
as the control-input signal. The signal choice satisfies the 
observability/controllability index based on modal analysis 
(e.g., residues analysis method [16]). 

e The important part of WADC is the dynamic output-feedback 
control matrix Kpor. Unlike the traditional control gain of 
output-feedback control, the Kpor is not a constant value but 
has a state-space (SS) or a transfer-function (TF) representa- 
tion. It can provide phase compensation for the concerned LFO 
mode as well as delay compensation for the wide area signal. 

e The high-pass filter (HPF) is used to preprocess the control signal 
and retain the concerned LFO mode that should be damped. The 
time-constant Twn selects 10 s according to the experience [17]. 


In general, the objective of HVDC DOF-WADC is to damp LFO 
especially the inter-area oscillations in power systems and at the 
same time ensure the required control performance of the HVDC 
local pole controller. 


3. Controller design considering time-varying delay of wide 
area signal 


3.1. Linearized discrete-time model with time-varying delay 


The linearized discrete-time model, which takes the time- 
varying delay of wide area signal into consideration, can be 
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Fig. 1. Basic structure of HVDC wide area damping control. 


represented as 


x(k+ 1) = Ax(k)+ Buk), 
y(k) = Cgx(k—d(k)), (1) 
x(k) = p(k), —d2<k<0O, 


where x(k) is the state vector; u(k) and y(k) are the wide area 
control-output and -input, respectively; A, B, and Cy are the 
system, the input and the output matrices, respectively; d(k) is 
the time-varying delay that satisfies 


dı < dk < d2. (2) 
The DOF controller can be represented as follows: 
X-(k+1) _ Ac Be | | Xc(k—d(k)) 3 
u(k) | Ce De y(k) : 6) 
where xk) is the state vector of DOF controller; Ae Be Ce and De 
are the real matrices. 


Under the action of (3), the open-loop system (1) yields a 
closed-loop system as follows: 


élk+1) = AE(k) + Bae(k—d(k)), (4) 
where 

xb] . TA O s BD:Ca BCe 
w= w| =[o of B= | ac, | 


The By in (4) can be further described as 
g [2 0]|Pe Ce|fCa 0]_fB Olp Ca 0 
=o | Bo Ac E |= [o ‘| sor| o ik 


The optimization of DOF controller in (3) is transformed into 
that of static output-feedback (SOF) controller Ksor, that is 


> De Ce 
u(k) = K sory(k) = Be Ac Wk). (5) 


The aim of controller design is to find a suitable Ksor to ensure 
the closed-loop system (4) can be asymptotically stable even 
when the control signal has a certain time-varying delay. 


3.2. Stabilization conditions 


Currently, there are several stabilization analysis methods 
generalized in LMI framework [18]. However, these methods 
cannot solve the effects of time-varying delay on stabilization 
performance effectively. In this paper, based on free-weighting 
matrix method [15] developed recently for time delay system, a 
theorem will be derived to analyze the stability of the system 
(1) with time-varying delay of wide area signal. 


Theorem. For given scalars d4 and dz with dz > d, > 0, system (1) 
with dp satisfying (2) can be stabilized by (3) if there exist P= P" > 0, 
L=L'>0, Q=Q/>0 (i=1, 2, 3) Z=Z)>0, Rj=R >0 
(G=1, 2) X=Řp X22 =e y2]>0,N=[N] N3, S=[S; S3], 
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M =[M] M}, and Řsor satisfies that 


1 oa (0) 6 
Ww) = | yr z,| 2% (6) 
2 o (0) 7 
Y )= sT Z 2U, ( ) 
3 X+Y M 6 3 
Y )= MT Zi +23 2U, ( ) 
PL=I, Z;Rj=I1G=1,2), (9) 
o> =Æ bZ dē 
51 —L 0 0 
0, 10 
d=, 0 -dR 0 |~ (0) 
4252 (0) 0 —dy2R2 


od n2 = Sr -Mı 
2 2 ¢ S% -M 
Ps Ss oy oS 
-Mi -M 0 -Q 


pi = Q14+Q2+(di2 + 1)Q3—P+Ni +N] +d2X11 +d12Y11, 


oy2= N35—N, +M,—S; + 42X12 +d12Y12, 


bon = —Q3—N2—N5 + Mp + Mj—S2—S5 + d2X22 +d12Y 22, 
Z1 =[A BKsorCa 0 0), 

E =|A—I BK sopCq 0 0), 

di2 = dy—d. 


Proof. Define 
nD = 1+ 1)—c(). (11) 


It is clear from (11) that €(k+1)=€(k)+n(k) and 
n(k) = E(k+1)—€(k). Then, choose a Lyapunov functional candi- 
date to be 

Vk) = Vi (k) + V2(k) + V3(k) + Va(k), 

Vi (k) = E*PE(K), 


k-1 k-1 
Vaky= So EaD J EW2EW, 
i=k-d i= k-dy 
kad) k-1 (12) 
v= X Sle QE), 
l= k-dzi =l 
—1 k—1 —dı—1 k-1 
Vaky= X do n@Zind+ X $ n OZO, 
i=-dl=k+i i=-—dl=k+i 


where P= P" >0,Q;=Q7 > 0 (i=1,2,3),Zi=Z7 > 0 (i=1,2). 
Defining AV(k) = V(k+1)—V(k), calculating the AV(k) along the 

system (1) and based on the free-weighting matrix method, 

following the similar analysis process in [15], it yields: 


CK b+ IPE + F5[doZ1 +d1222122} 6) 
= Thay C2 (KD Pi ék D- eg! GD PED 
sy 2D P3E2(kD 


AV(k) < 


where ¢(t) = [E (OE t-d) é" (t-d1)é" (t-d), Gr (kK) = [E OE t- 
Tk), Colk) = [éi don. 

Therefore, if Y; >0(i= 1,2,3) and ġ+21P21 +2}3(d2Z1 +d1222) 
E2 > 0, then AV(k) <0. By using Schur complements, it is found 
that p+ 2) PE, + 25(doZ, +d12Z2)Z> > 0 is equal to 


o ETP do55Z, dZ 

P=, -P 0 0 
ie, o0 -daz o |=" a3) 
dy2Z55, 0 0 di2Z2 


For the inequality (13), pre- and post-multiplying its left and 
right sides by diag{1,P-',Z7',Z3' } yields 


$ E dZ dē} 

=, —P! 0 0 
Ga: 0 stg? o |<® ue) 
d2% 0 0 dy2Z5' 


Defining L=P~!, R;=Z;'(i=1,2) in (14), it can be seen that 
(9) and (10) hold. Therefore, based on the above analysis, it can be 
concluded that if conditions (6)-(10) are satisfied, system (1) can 
be stabilized by the controller (3). 


3.3. Optimization objective and iterative algorithm 


Since there are some equality terms like (9) in the theorem, it 
is impossible to optimize the controller (5). To overcome this 
problem, the cone complementary linearization (CCL) algorithm 
[19] is used to transform the controller design into a nonlinear 
optimization problem as follows: 


2 
Min. Tr (re 5 zn) 


j=1 


(6)—(8),(10),and 


i p| 20 Z I| o, j=1,2 
stl L PJS? JE Re (15) 


P>0, L>0, Q;>0, i=1,2,3, 
Z; > 0, Rj >0, j=1,2. 


The following iterative algorithm is planned to solve the above 
nonlinear optimization problem. 


Algorithm. 


è Step-1: Set the initial value dz > d4, the iteration step Ad and 
the iteration number N; 

e Step-2: Find a set of initial values (Po, Lo, Qio, Zjo, Rjo, Mo, No, So, 
Xo, Yo, Ksoro) so that (15) can be held; 
Set k=0; 

e Step-3: Solve the following optimization problem: 


2 
Min.Tr (aa + So ZiRit sza) 
j=1 


s.t. (15) 
Set d$ =d; 

e Step-4: Make 
Rik+1 =R; G=1,2). 


Prat =P, Lk+1 =L Zjk+1 =Z;, and 
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Input A, B, Ca, 
Set initial delay dz > dj, iteration step Ad 
and iteration number N. 


Find a set of initial values (Fo, Lo, Qio, Zjo, Ryo, Mo, 
No, So, Xo, Yo, Ksoro) so that (15) can be held; 
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Set k = 0. 


s. t. (15) 
Set ds = a. 


Min. Tr (PL+ PLs + Y3 (Zj Rix + ZjeRs)) 


y 


Make Pk+1 = Po, Lk+1 = Lo, Zijk+1 = Zj and Rj,k+1 = Rj, (j = 1,2). 


Satisfy (6)-(8) and 
$ 


=&PT -P 
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Output Ksor 


Fig. 2. Flowchart of the optimization algorithm for controller design. 


è Step-5: If (6)-(8) can be held and the following LMI is feasible: 
$ ZIP do55Z, di2E5Z2 


5,P’ -P 0 0 ; 
dZ 0 -dZ o | 
d1222Z} 0 0 dy2Z2 


then, set d5+! = dk +Ad; go to Step-2; 
else if k >N; 
stop the iteration; 
else 
k=k+1; go to Step-3; 
end 
end 
e Step-6: Output the final results; End. 


Fig. 2 shows the flowchart of the controller optimization. From 
this figure, it is clear that the time-varying delay of the wide area 
signal is considered during the controller design, which is helpful to 
find a suitable controller that can keep on an effective damping under 
the condition of time-varying delay on the wide area control signal. 


4. Case-I: four-machine two-area test system with HVDC/AC 
interconnected lines 


The four-machine two-area test system [17] is used to validate 
the presented controller design method. This system is divided 


DOF-WADC |«------- 


Fig. 3. Four-machine two-area test system with HVDC/AC interconnected lines. 


into two areas where G1 and G2 are located in one area and G3 
and G4 are located in the other area, as shown in Fig. 3. There are 
parallel AC/DC transmission lines to interconnect these two areas. 

This test system has a critical inter-area oscillation mode, 
which is represented as the generator oscillations between G1-G2 
and G3-G4. Here, to damp this mode, the supplementary control 
of HVDC transmission is utilized to construct the DOF-WADC 
strategy. The control-input chooses the power flow in lines 9-10 
and the control-output is the supplementary control-input of the 
converter pole controller at the rectifier side of HVDC transmis- 
sion, as shown in Figs. 1 and 3. 

A line-to-ground fault near bus 9 of lines 8-9, which starts at 
1.0s and is cleared after 0.1 s, is simulated to illustrate the 


5752 Y. Li et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5747-5759 


damping performance of the designed HVDC DOF-WADC when 
system suffers this big disturbance. Fig. 4 shows the dynamic 
response of the rotor speeds of the generators located in different 
areas. It can be seen that for the system without HVDC DOF- 
WADC, a big disturbance excites serious power oscillations 
among generators, as shown in Fig. 4(a). But when the system is 
installed with the HVDC DOF-WADC, the oscillations are damped 
effectively, as shown in Fig. 4(b). 

Furthermore, Fig. 5 represents the dynamic response of the 
relative angle between G1 located in one area and G3 located in 
the other area. Fig. 6 represents the control-output response of 
the HVDC DOF-WADC. From these, it is clear that by means of the 
HVDC DOF-WADC, the inter-area oscillations are damped greatly. 
Besides, DOF-WADC can quickly respond to the oscillations and 
provide the effective control signal for oscillations damping. 

Fig. 7 shows the dynamic response of power flow on tie-lines. 
From this, it can be seen that the big disturbance leads to serious 
power oscillations on the backbone transmission lines, which 


a 1.025 


1.02 
1.015 
1.01 
1.005 


Generator speed/pu 


oO 


Generator speed/pu 


0.9995 


5 10 


further makes system collapse. But when using HVDC DOF- 
WADC, the power oscillations are damped effectively. Therefore, 
the HVDC DOF-WADC strategy is good to the stable and security 
operation of interconnected system. 

To investigate the effects of time-varying delay of wide area 
signal on the damping performance, various time delays, includ- 
ing the transmission delay and the signal-processing delay, are 
considered as a total delay. Fig. 8 shows the damping perfor- 
mance of the HVDC DOF-WADC under different time delays. It is 
clear that the designed controller can maintain good damping 
performance, though there exists a certain time-varying delay. 


5. Case-II: sixteen-machine five-area test system with HVDC/ 
AC interconnected lines 


To further validate the effectiveness of the presented DOF- 
WADC, a more complex sixteen-machine five-area test system, 


15 20 25 30 


Time/s 


Fig. 4. Dynamic response of generator speed: (a) without HVDC DOF-WADC; (b) with HVDC DOF-WADC. 
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Fig. 5. Dynamic response of the relative angle between G1 and G3. 
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Fig. 6. Control-output of HVDC DOF-WADC. 
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Fig. 7. Dynamic response of power flow: (a) lines 10-11; (b) lines 6-7. 


that is the New England Test System (NETS) and the New York 
Power System (NYPS) interconnected system [17], shown in Fig. 9, 
is studied in detail. 


5.1. Choice of optimal wide area signal for DOF-WADC 


At the first, the modal analysis is carried out to reveal the 
inter-area oscillations. Fig. 10 shows the oscillation shape and the 
eigenvalues. On the nominal operating condition, this system has 
four typical inter-area modes, as shown in Fig. 10(b). Among 
these, there is a mode 1 with frequency f = 0.5773 Hz and small 
damping ratio p = 0.0152. Fig. 10(a) shows its oscillation shape, 
from which it can be seen that this mode is dominated by the 
oscillations between area 3 and areas 1 and 5. 

The HVDC DOF-WADC is used to damp this critical mode. Here, 
power flow is candidate as the control-input. According to residue 
method, the damping contribution of power flow on each line of 


the test system can be obtained, as shown in Fig. 11. It can be see 
from this figure that compared to the local control signals, the 
wide area signal, that is the power flow in line 85 (see lines 52-68 
in Fig. 9), can get the best control performance. Therefore, it is 
selected as the control-input of DOF-WADC. 


5.2. DOF-WADC design 


The order of the linearized model of this test system is 202, 
which is too high to be difficult for controller design within the 
LMI framework. Here, the balanced truncation model-reduction 
method is used to reduce the model order from 202 to 9. Fig. 12 
shows the bode diagram of the full- and the reduced-order system. 
It is clear that over the wide frequency range (10-?—10* Hz), the 
reduced-order system is very close to the full-order system, which 
indicates that the reduced-order system can fully reflect the 
dynamic characteristics of the full-order system. 
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Fig. 8. Delay effect of wide area signal on damping performance of HVDC DOF-WADC. 
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Fig. 9. Sixteen-machine five-area NETS-NYPS interconnected system with HVDC/AC interconnected lines. 


According to the presented controller design method, 
the HVDC DOF-WADC for this test system can be obtained 


as follows: 


De 
u(k) = lz 


Ce i 
Ac yk) = | 


0.0282 0.6303 
—0.0024 0.9862 


|y% 


(16) 


Furthermore, Fig. 13 gives the root locus of the test system 
which adopts the above controller (16) with different overall 
control gain. From this figure, it can be seen that the designed 
controller represents a good response performance to the con- 
sidered inter-area mode. It can greatly move this mode (small 
damping ratio, p <0.05) to the left-plane, which means the 
improvement of the system stability. 
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Fig. 10. Inter-area mode with small damping ratio (p < 0.05): (a) oscillation shape; (b) eigenvalues. 
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Fig. 11. Damping contribution of different line flow as the control-input of DOF-WADC. 
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Fig. 12. Comparison of bode diagram of the full- and the reduced-order system. 
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5.3. Nonlinear simulation 


In order to evaluate the damping performance of the presented 
DOF-WADC, the nonlinear simulation is performed, which 
includes the controller evaluation on different operating condi- 
tions (e.g., line-to-ground fault and line outage). 


5.3.1. Line-to-ground fault 
A line-to-ground fault near bus 1 of lines 1-2 is considered to 
test the performance of the designed HVDC DOF-WADC. Fig. 14 
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shows the dynamic response of the relative angle between gen- 
erators located in different areas. From this figure it can be seen 
that the designed controller can effectively damp the inter-area 
oscillations within 15 s. Fig. 15 gives the control-output, which also 
indicates that the controller has a quick response for oscillation 
damping. Furthermore, Fig. 16 represents the power flow on the 
tie-lines, which indicates that the power oscillations on the back- 
bone transmission lines are damped effectively. 

Besides, to evaluate the damping performance under time- 
varying delay of wide area signal, Fig. 17 further gives the 
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Fig. 13. Root locus diagram of the test system under the action of the DOF-WADC. 


Relative angle between G15 and G16/° 


oO 


Relative angle between G14 and G16/° 


Without DOF-WADC 
With DOF-WADC 


* Without DOF-WADC 
With DOF-WADC 


15 20 25 30 
Time/s 


Fig. 14. Dynamic response of relative angle: (a) between G15 and G16; (b) between G14 and G16. 
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Fig. 15. Control-output response of DOF-WADC to a line-to-ground fault. 
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Fig. 16. Dynamic response of power flow in the tie-lines: (a) lines 41-66; (2) lines 52-68. 
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Fig. 17. Time-varying delay effect of wide area signals on damping performance of DOF-WADC. 
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Fig. 18. Dynamic response of the test system to a line outage: (a) power flow in lines 41-66; (b) relative angle between generators. 
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Fig. 19. Dynamic response of the test system under different time delays of wide area signal. 


dynamic response of the power flow in lines 52-68. It is clear 
from this figure that although it has different time delays for the 
wide area signals, the designed DOF-WADC still maintains a good 
damping performance, which illustrates the effectiveness and the 
reliability of the presented controller design method. 


5.3.2. Line outage 

The outage of interconnected lines 1-47 is considered to 
further test the DOF-WADC under different operating conditions. 
Fig. 18 shows the dynamic responses of the test system to such 
line outage. From both the power flow from the transmission side 
and the relative angle from the generation side, it is clear that the 
designed controller still keeps on a good damping performance 
under this operating condition. Moreover, Fig. 19 gives the 
dynamic responses of the test system under different time delays 
of wide area signals. It is also clear that the DOF-WADC maintains 


a good damping performance even when the wide area signal 
exists various time delays. 


6. Conclusion 


This paper proposed a DOF-WADC and related controller 
design method to damp inter-area oscillations and thus to 
enhance the stability of large-scale power systems with parallel 
HVDC/AC interconnected lines. This method can sufficiently 
utilize the supplementary control function of HVDC transmission, 
and implement an effective wide area damping control through 
introducing a suitable wide area signal. By employing an 
advanced stabilization theorem, the design problem of DOF- 
WADC is formulated as a standard nonlinear optimization pro- 
blem with a set of LMI constrains. An iterative algorithm is 
correspondingly presented to solve the controller parameters 
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within the LMI framework. The time-varying delay of wide area 
signals is considered during the process of controller design. 
Unlike the conventional SOF control that only has a single 
constant control gain, in this paper, the WADC has a good 
dynamic damping performance, and it is represented as a 
dynamic SS form. Two case studies validated this control method, 
and revealed the good damping performance when applied in 
power systems with parallel HVDC/AC interconnected lines. 
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